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Abstract 

We present a theoretical framework that accounts for the new Dj and D s j mesons measured 
in the open-charm sector. These resonances are properly described if considered as a mixture of 
conventional P— wave quark-antiquark states and four-quark components. The narrowest states 
are basically P— wave quark-antiquark mesons, while the dominantly four-quark states are shifted 
above the corresponding two-meson threshold, being broad resonances. We study the electromag- 
netic decay widths as basic tools to scrutiny their nature. The proposed explanation incorporates 
in a natural way the most recently discovered mesons in charmonium spectroscopy. 
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I. INTRODUCTION 



During the last few years, heavy meson spectroscopy is living a continuous excitation 
due to the discovery of several new charmed mesons. Two years ago BABAR Collaboration 
at the Standford Linear Accelerator Center (SLAC) reported the observation of a charm- 
strange state, the _D*j(2317) 0. It was confirmed by CLEO Collaboration at the Cornell 
Electron Storage Ring and also by Belle Collaboration at KEK @. Besides, BABAR 
had also pointed out to the existence of another charm-strange meson, the -D^j (2460) hi. 
This resonance was measured by CLEO and confirmed by Belle Belle results p[ 
are consistent with the spin-parity assignments of J p = + for the Z)*j(2317) and J p = 
1 + for the D s j(2460). Thus, these two states are definitively well established, confirmed 
independently by different experiments. They present unexpected properties, quite different 
from those predicted by quark potential models. If they would correspond to standard 
P— wave mesons made of a charm quark, c, and a strange antiquark, s, their masses would 
be larger Q, around 2.48 GeV for the D* sJ (2317) and 2.55 GeV for the D aJ (2A60). They 
would be therefore above the DK and D*K thresholds, respectively, being broad resonances. 
However the states observed by BABAR and CLEO are very narrow, F < 4.6 MeV for the 
£^(2317) and F < 5.5 MeV for the L> sJ (2460). 

The intriguing situation of the charm-strange mesons has been translated to the non- 
strange sector with the Belle observation of a nonstrange broad scalar resonance, Dq, 
with a mass of 2308 ± 17 ± 15 ± 28 MeV/c 2 and a width r = 276 ±21 ±18 ±60 MeV. A state 
with similar properties has been suggested by FOCUS Collaboration at Fermilab during 
the measurement of masses and widths of excited charm mesons D^ ■ They found significant 
evidence for a broad excess parametrized as an S— wave resonance. Without being able to 
clearly distinguish the origin of this broad excess, they conclude that their results are in 
agreement with the Belle results [Hf. This state generates for the open-charm nonstrange 
mesons a very similar problem to the one arising in the strange sector with the _D*j (23 17). 
If the /}q(2308) would correspond to a standard P— wave meson made of a charm quark, c, 
and a light antiquark, n, its mass would have to be larger, around 2.46 GeV. In this case, 
the quark potential models prediction and the measured resonance are both above the Dtt 
threshold, the large width observed being expected although not its low mass. 

The last step in this series of discoveries has been the observation of a new charm-strange 
meson, D s j, announced by SELEX Collaboration at Fermilab p] with a mass of 2632.5 ±1.7 
MeV/c 2 and a small width, T < 17 MeV. However, up to now no other experiment has been 
able to confirm the existence of this resonance 0. 

There have been many theoretical interpretations for the masses and widths of the new 
resonances, but most part of them have been devoted to explain the strange states. Ref. 
made use of a unitarized meson model, the existence of a quasi-bound cs state due to 
the coupling with the nearby S— wave DK threshold. The coupling to the DK channel 
within a QCD string model an a chiral Lagrangian was used in Ref. [Hi • Ref. [IH proposed 
a cs structure with modified non-central forces. Ref. combined HQET with a chiral 
effective Lagrangian to interpret these states as the missing j = 1/2 member of the cs L = 1 
ground state multiplet, being j the angular momentum of the strange quark. The smaller 
mass of the _D*j (23 17) is attributed in Ref. [lj| to a coupled channel effect. Finally more 
involved qualitative solutions like DK— molecules 3|> D s n— atoms jl^], four-quark states 



16l 1 171 ] and their combination with qq states |18[ have been invoked to explain why the 



new states have a nature so different from canonical qq mesons. Although the nonstrange 
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partners of the D s /s have received much less attention, as discussed above, they present 
similar spectroscopic properties that should be acknowledge altogether with those of the 
strange states in any reliable model. 

The difficulties to identify the Dj and D sJ states with conventional cq mesons are rather 
similar to those appearing in the light-scalar meson sector and may be indicating that 
other configurations are playing a role, as could be for example four-quark contributions. 
qq states are more easily identified with physical hadrons when virtual quark loops are 
not important. This is the case of the pseudoscalar and vector mesons, mainly due to the 
P— wave nature of this hadronic dressing. On the contrary, in the scalar sector it is the qq 
pair the one in a P— wave state, whereas quark loops may be in a S— wave. In this case 
the intermediate hadronic states that are created may play a crucial role in the composition 
of the resonance, in other words unquenching may be important. The vicinity of these 
components to the lightest qq state implies that they have to be considered. This has been 
shown as a possible interpretation of the low-lying light-scalar mesons, where the coupling 
of the scalar qq nonet to the lightest qqqq configurations allows for an almost one-to-one 
correspondence between theoretical states and experiment [20|. 

In this work we pretend to explore the same ideas for the understanding of the properties 
of the Dj and D s j meson states. For this purpose, in the next section we will present our 
calculating scheme and the basic ingredients of the constituent quark model used. Section 
III will be devoted to present and discuss our results in connection with those obtained for 
the light-scalar mesons. Finally, in Sec. IV we will resume the most important conclusions 
of our work. 

II. CALCULATING FRAMEWORK 

In non-relativistic quark models gluon degrees of freedom are frozen and therefore the 
wave function of a zero baryon number (B=0) hadron may be written as 

\B = Q) = tt l \qq) + tl 2 \qqqq) + .... (1) 

where q stands for quark degrees of freedom and the coefficients Qi take into account the 
mixing of four-quark and qq states. |B = 0) systems could then be described in terms of a 
hamiltonian 

H = H + H 1 being H =( H « ) H t = ( ° , ( 2 ) 

V U El qqqq J \ ' qq^qqqq U / 

where Hq is a constituent quark model hamiltonian described below and Hi, that takes 
into account the mixing between qq and qqqq configurations, includes the annihilation op- 
erator of a quark- ant iquark pair into the vacuum. This operator could be described using 
the 3 Po model, however, since this model depends on the vertex parameter, we prefer in 
a first approximation to parametrize this coefficient by looking to the quark pair that is 
annihilated and not to the spectator quarks that will form the final qq state. Therefore we 
have taken V qq ^ qqqq = 7. If this coupling is weak enough one can solve independently the 
eigenproblem for the hamiltonians H qq and H qqqq , treating Hi perturbatively. To ensure that 
the perturbative treatment is justified, 7 cannot take all possible values, being restricted to 
\j/(E^ PC - Ej£c)\ 2 < 1. This restriction will limit the energy range of the mixed states 
once the unmixed energies are calculated. The two-body problem has been solved exactly by 



3 



means of the Numerov algorithm |21|. The four-body problem has been solved by means of 
a variational method using the most general combination of gaussians as trial wave functions 
[22! . In particular, the so-called mixed terms (mixing the various Jacobi coordinates) 
that are known to have a great influence in the light quark case have been considered. 

Although the constituent quark model used is described in Ref. let us outline here 
its basic ingredients. Since the origin of the quark model hadrons have been considered to be 
built by constituent (massive) quarks. Nowadays it is widely recognized that the constituent 
quark mass of light quarks appears because of the spontaneous breaking of the original chiral 
symmetry of the QCD Lagrangian, what gives rise to boson-exchange interactions between 
quarks. The different terms of the potential contain central and tensor or central and spin- 
orbit contributions that will be grouped for consistency. Therefore, the chiral part of the 
quark-quark interaction can be resumed as follows, 



v x(ra) + V?(nj) + v x so (f tJ ) , (3) 

where C stands for central, T for tensor, and SO for spin-orbit potentials. The central part 
of the quark-quark meson-exchange potentials are given by: 
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the angle 9p appears as a consequence of considering the physical 77 instead the octet one. 
g c h = m q /f n , the X's are the £77(3) flavor Gell-Mann matrices and the cr's are the spin 
quark Pauli matrices, is the quark mass and m n , m^ and m v are the masses of the £77(3) 
Goldstone bosons, taken to be their experimental values. The Aj's are cutoff parameters. m a 
is determined through the PCAC relation m 2 a ~ m\ + Am 2 ud . Finally, Y(x) is the standard 
Yukawa function defined by Y(x) = e~ x /x. 

There are three different contributions to the tensor potential, 
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being SV,- = 3 [Bi • fij){Bj ■ r^) — o% ■ Bj the quark tensor operator and H{x) = (1 + 3/x + 
3/x 2 )Y{x). 

Finally, the spin-orbit potential only presents a contribution coming form the scalar part 
of the interaction, 
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where G(x) = (1 + 1/x) Y(x)/x. 

QCD perturbative effects are taken into account through the one-gluon-exchange (OGE) 
potential The nonrelativistic reduction of the one-gluon-exchange diagram in QCD for 
point-like quarks presents a contact term that, when not treated perturbatively, leads to 
collapse [25]. This is why one maintains the structure of the OGE, but the 5 function is 
regularized in a suitable way. This repolarization, justified by the finite size of the systems 



studied, has to be flavor dependent |2( 
reads, 
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where A c are the SU(3) color matrices, a s is the quark-gluon coupling constant, and r (/i) = 
rol^nn/Hij, where //^ is the reduced mass of quarks ij (n stands for the light u and d quarks) 
and f is a parameter to be determined from the data. 

The noncentral terms of the OGE present a similar problem. For point-like quarks they 
contain an 1/r 3 term. Once again the finite size of the constituent quarks allows for a 
regularization, obtaining tensor and spin-orbit potentials of the form, 
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where S± — Si ± Sj, and r s (/x) = r g j\i presents a similar behavior to the scaling of the 
central term. 

The strong coupling constant, taken to be constant for each flavor sector, has to be scale- 
dependent when describing different flavor sectors j23|. Such an effective scale dependence 
has been related to the typical momentum scale of each flavor sector assimilated to the 
reduced mass of the system 28||. This has been found to be relevant for the study of the 
meson spectra and parametrized in Ref. (2lJ as 
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where \i is the reduced mass of the interacting qq pair and ao, Ha, and 70 are fitted param- 
eters. 

Finally, any model imitating QCD should incorporate confinement. Lattice calculations 
in the quenched approximation derived, for heavy quarks, a confining interaction linearly 
dependent on the interquark distance. The consideration of sea quarks apart from valence 
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quarks (unquenched approximation) suggests a screening effect on the potential when in- 
creasing the interquark distance j2^|. Creation of light-quark pairs out of vacuum in between 
the quarks becomes energetically preferable resulting in a complete screening of quark color 
charges at large distances. String breaking has been definitively confirmed through lattice 
calculations [30] in coincidence with the quite rapid crossover from a linear rising to a flat 
potential well established in SU(2) Yang-Mills theories [3j|. The central part of a screened 
potential simulating these results can be written as, 

Vg ON {nj) = -oe (1 - e~^)(X\ • &j) ■ (12) 

At short distances it presents a linear behavior with an effective confinement strength a = 
a c fi c X c i ■ X c j, while it becomes constant at large distances. Screened confining potentials 
have been analyzed in the literature providing an explanation to the missing state problem 
in the baryon spectra j^l, improving the description of the heavy-meson spectra [33|, and 
justifying the deviation of the meson Regge trajectories from the linear behavior for higher 
angular momentum states [34| . 

Confinement also presents an spin-orbit contribution taken to be an arbitrary combination 
of scalar and vector terms of the form 

VS8 N (r,) = -CK ■ Al) [(K 2 + m)){l - 2 a.) 

+4m im ,(l - a,)) (S + ■ L) + (m) - m?)(l - 2a s )(S_ • L)] (13) 

where a s would control the ratio between them. 

Once perturbative (one-gluon exchange) and nonperturbative (confinement and chiral 
symmetry breaking) aspects of QCD have been considered, one ends up with a quark-quark 
interaction of the form 

t / / QiQj = nn/ sn V CO n + Voce + V x , . 

mj \q l q ] = cn/cs ^V CON + Voce ' 1 } 

Note that for the particular case of heavy quarks, chiral symmetry is explicitly broken and 
therefore Goldstone boson exchanges do not contribute. The model parameters and a more 



detailed discussion of the model can be found in Refs. 120, 21, 2 



III. RESULTS AND DISCUSSION 



A thoroughly study of the full meson spectra has been presented in Ref. |21(. The results 
for the open-charm mesons are resumed in Table |1J It can be seen how the open-charm 
states are easily identified with standard cq mesons except for the cases of the _D*j(2317), 
the _D s j (2460), and the Dq(2308). We will also comment on the recent measurement of 
the -D l j ) (2430). This behavior is shared by almost all quark potential model calculations jij. 
Although the situation from lattice QCD is far from being definitively established, similar 
problems are observed. Lattice NRQCD in the quenched approximation predicts for the 
_D*j(2317) a mass of 2.44 GeV j36[, while using relativistic charm quarks the mass obtained 
is 2.47 GeV [3711 . U nquenched lattice QCD calculations of cs states do not find a window for 
the _D*j(2317) [19(, supporting the difficulty of a P— wave cs interpretation. The quenched 
lattice QCD calculation of the spectrum of orbitally excited D s mesons of Ref. [38| concludes 



6 



that, although the results obtained are consistent with a cs configuration, the statistical and 
systematical uncertainties are too large to exclude the exotic states based on potential quark 
models |4| . The same situation may be drawn from heavy quark symmetry arguments. One 
finds that the scalar cs state belongs to the j — 1/2 doublet, but since the j = 3/2 doublet 
is identified with the narrow _D S 2(2573) and D s i(2536) (with total widths of 15^4 MeV and 
< 2.3 MeV, respectively) the scalar state is expected to have a much larger width than the 
one measured for the _D*j(2317) [3ot | . 

Thus, one could be tempted to interpret these states as four-quark resonances within the 
quark model. The results obtained for the cnsn and cnnn configurations with the same 
interacting potential of Sec. II are shown in Table |HJ The 1=1 and 1 = states are far 
above the corresponding strong decay threshold and therefore should be broad, what rules 
out a pure four-quark interpretation of the new open-charm mesons. 

As outlined above, for P— wave mesons the hadronic dressing is in a S— wave, thus physical 
states may correspond to a mixing of two- and four-body configurations, Eq. (JTJ). In the 
isoscalar sector, the cnsn and cs states get mixed, as it happens with cnnn and cn for the 
1 = 1/2 case. Let us notice that the interacting potential contains terms mixing P— wave 



S = and S = 1 states [21[, arising from the different spin-orbit contributions for systems 
made of quarks of different mass. The mixing angle comes determined by the structure of 
the interaction. The parameter 7 has been fixed to reproduce the mass of the D*j(2317) 
meson, being 7 = 240 MeV. Using this value one has \^/(E^ PC - Ej^c)\ 2 ~ 0.25 < 1, 
a ratio consistent with the assumption that the hamiltonian H\, Eq. (J2J), can be treated 
perturbatively. This perturbative condition together with \j/E gq -\ « \j/E qq gq\ w 0.1 < 1 
also ensures that the mixing with higher states can be neglected. 

The results obtained are shown in Table ITTT1 Let us first analyze the nonstrange sector. 
The 3 Po cn pair and the cnnn have a mass of 2465 MeV and 2505 MeV, respectively. Once the 
mixing is considered one obtains a state at 2241 MeV with 46% of four-quark component 
and 53% of cn pair. The lowest state, representing the -Dq(!2308), j s above the isospin 
preserving threshold Dir, being broad as observed experimentally. The mixed configuration 
compares much better with the experimental data than the pure cn state. The orthogonal 
state appears higher in energy, at 2713 MeV, with and important four-quark component. 
A similar process would modify the state representing the _D ( j ) (2430), but in this case one 
would need the mass of the / = 1/2 J p = 1 + four-quark state. The huge basis generated for 
such quantum numbers makes unfeasible its calculation. A correct description of this state 
would require a mass of 2.9 GeV for the four quark state mentioned above. 

Concerning the strange sector, the _D*j(2317) and the _D s j(2460) are dominantly cs J = 
+ and J = 1 + states, respectively, with almost 30% of four-quark component. Without 
being dominant, it is fundamental to shift the mass of the unmixed states to the experimental 
values below the DK and D*K thresholds. Being both states below their isospin-preserving 
two-meson threshold, the only allowed strong decays to D*tt would violate isospin and are 
expected to have small widths. This width has been estimated assuming either a qq structure 
T2I 4oj . a four-quark state 41 1 or vector meson dominance obtaining in all cases a width 



of the order of 10 keV. The second isoscalar J p = 1 + state, with an energy of 2555 MeV 
and 98% of cs component, corresponds to the -D s i(2536). Regarding the _D*j(2317), it has 
been argued that a possible DK molecule would be preferred with respect to an I = cnsn 
tetraquark, what would anticipate an / = 1 cnsn partner nearby in mass |l4j|. Our results 
confirm the last argument, the vicinity of the isoscalar and isovector tetraquarks (see Table 
in}, however, the restricted coupling to the cs system allowed only for the 1 = four-quark 
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states opens the possibility of a mixed nature for the D* t/ (2317), the 1 = 1 tetraquark 
partner remaining much higher in energy. The I = 1 J = + and J = 1 + four-quark states 
appear above 2700 MeV and cannot be shifted to lower energies. 

Our results do not show any state around 2600 MeV, the mass region compatible with 
the _Dj}(2632) measured by SELEX. Its two possible theoretical partners, the 2 3 S' 1 state 
(suggested in Ref. 0) or the 1 3 D 1 state, lie above 2700 MeV (2764 and 2873 MeV, respec- 
tively). This state has been proposed as the csss partner of the _D*j(2317) [T3|. Naively, 
the _D*j (23 17) sets a scale of 2320 MeV for the cnsn sector, and an augment of 150 MeV 
for each strange quark could accommodate a csss system near the mass of the _Dj}(2632). 
However, the mass of the _D*j(2317) is obtained through the coupling to the cs quark pair, 
being the mass of the tetraquark configuration above 2700 MeV, disregarding a possible 
csss interpretation of the _Dj}(2632). It has also been discarded as the first radial excita- 
tion of the -D *(21 12) A careful analysis of several theoretical interpretations has been 
done in Ref. |44j, the surprising properties of this state not fitting in any of the scenarios 
considered. Several experiments have failed trying to confirm the existence of this state . 
The confirmation or refutation of the _Dj}(2632) is clearly an important priority for meson 
spectroscopy. 

The structure of the _D*j(2317) and the _D s j(2460) mesons could be scrutiny, apart from 
their masses, also through the study of their electromagnetic decay widths. Using the 
standard formalism described, for example, in Ref. jiol ]. the El radiative transitions for 
cs — > cs + 7 are given by: 
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r[ cs - cs + 7 ] = -a (c Q )' z uj\2J f + 1) ( 2S+1 Sj, \r\ 2S+1 Pj) ' S tf (15) 

where Sif is a statistical factor with Sif = 1 for the transitions between spin-triplet states 
and Sif = 3 for the transition between spin-singlet states, (eg) is an effective quark charge 
given by 

eg = "\ 16 

where Ci is the charge of the quark(antiquark) i in units of |e|, a is the fine structure constant 
and io is the photon energy. Once the mixing between two- and four-quark components 
has been included there would also be a contribution to the decay width coming from 
cnsn — > cs + 7, that goes necessarily through the annihilation of a color singlet nn pair with 
photon quantum numbers within the four-quark wave function. The contribution to the 
electromagnetic decay width arising from the four-quark component gets suppressed. This 
can be illustrated by the analysis of the Ml electromagnetic decay width for cnsn(l + ) — > 
cnsn(0 + ) + 7, that is given by 

T[cnsn(l + ) -> cnsn(0 + ) + 7] = ~au 3 {2J f + 1) ((^o \fi z \ 3 S^ 2 (17) 

where 

4 ei 

From the above equation one obtains T[cnsn(l + ) — > cnsn(0 + ) + 7] m 4 eV, that will con- 
tribute to the r[D s j(2460) — > D*j(2317) +7] decay width, once the configuration probabil- 
ities given in Table ITTTl are considered, in the order of 0.25 eV. The value obtained in Ref. 
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[ioj for this decay, assuming a pure cs configuration, is one order of magnitude larger con- 
firming the smallness of the four to four-quark contribution. In Ref. |45j this decay width 
was estimated assuming a molecular structure, obtaining a larger value of the order of 17 
keV. This makes evident the important difference between the assumed molecular structure 
of Ref. 45j and the full multiquark configuration considered in this work. In the case of 
the four-quark configuration the orthogonality between the several components of the wave 
functions diminishes the decay widths, an effect that does not seem to be present within the 
molecular interpretation. The same reasoning applies to the T[cnsn — > cs + 7] decay width, 
being suppressed due to the reduced singlet- singlet component with photon quantum num- 
bers within the four-quark wave function. Similar conclusions were obtained in Ref. [l8| . 
As a consequence, the presence of a four-quark component diminishes the decay widths, 
making them different from those predicted for a pure cs state [3, 0, EHf ■ 

We compare in Table II VI our results for the radiative transitions of the _D*j (23 17) 
and _D s j (2460) with different theoretical approaches and the experimental limits re- 
ported by CLEO and Belle. The main difference is noticed in the suppression predicted 
for the -D sJ (2460) — > -D* + 7 decay as compared to the D sJ (2A60) -» D+7. A ratio 
.Dsj (2460) — > Dtl / Ds.t (24:60) — > D* +r y ^1 — 2 has been obtained assuming a qq struc- 
ture for both states [12, 40] (what seems incompatible with their properties). Heavy-hadron 
chiral perturbation theory calculations find a similar ratio 46]. We find a larger value, 
D sJ (2460) -> D+7/D sJ (2460) -> « 100, due to the small l 3 Pi cs probability of the 

D s j(2460). A similar enhancement has been obtained in Ref. ji^J (see penultimate column 
of Table HVj) in the framework of light-cone QCD sum rules in contrast to a previous calcula- 
tion of the same authors using vector meson dominance 0). As a consequence, the radiative 
transitions are an important diagnostic tool to understand the nature of these states. In 
view of the different predictions of the electromagnetic decay widths, a precise measurement 
of this decay would allow to distinguish not only between qq and non— qq states, but also 
between pure molecular and four-quark interpretations. 

Let us finally mention that the difficulties encountered for the interpretation of the new 
open-charm states as two-quark systems do not appear for the case of the recent charmed 
and B c states measured at different facilities. They do nicely fit into the predictions of the 
model used, see Table |VJ giving confidence to the results obtained in the present work. 

The interpretation we have just presented of the positive parity open-charmed mesons as 
a mixture of two- and four-quark states has also been used to account for the experimentally 
observed light-scalar mesons within the same constituent quark model 20], what gives us 
confidence on the mechanism proposed. Nonetheless, one should not forget that in the 
literature there is a wide variety of interpretations for the open-charm and also for the 
light-scalar mesons. As previously mentioned, they used to deal with one of both problems 
and even then with a particular set of states, being the most common ones the strange 
open-charm states and the isoscalar light mesons. 

Regarding the light-scalar mesons, there have been reported by the PDG two isovectors, 
five isoscalar and three isodoublet states. The situation is far from being definitively settled, 
neither from the experimental nor from the theoretical point of view. Experimentally, let us 
only mention that the recent analysis by BES of the J/ty — > <j)7i + 7i^ and J/^f — > <fiK + K~ 
data [53| requires a state /o(1790) distinct from the /o(1710). Recent reanalysis of the Crystal 
Barrel data 0| suggest the existence of a new state called /o(1200 — 1600). Theoretically, 
one would expect non— qq scalar objects in the mass range below 2 GeV. Multiquarks have 
been justified to coexist with qq states in the energy region around 1 GeV because they can 
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couple to ++ without orbital excitation |55|. Lattice QCD in the quenched approximation 
predicts the existence of a scalar glueball with a mass around 1.6 GeV 56l . 

From this complicated scenario many different interpretations of the light-scalar mesons 
have arisen. The overpopulation of ++ states below 2 GeV gave rise long ago to the 
speculation of the existence of four light-quark states. The most complete analysis was 
performed by Weinstein and Isgur 2^, 57], concluding that, normally, qqqq bound states do 
not exist, being the only exception the scalar sector where weakly bound states with a meson- 
meson molecule structure were found. Particular sets of states have been studied in the 
literature. There are several models analyzing the mixing between different configurations 
to yield the physical /o(1370), /o(1500), and /o(1710). Among them, Ref. [2p| assigning the 
larger glueball component to the /o(1710) is on the line with Refs. |58( and differ from those 
of Refs. jH^, 6(| concludin g th at the / (1710) is dominantly qq. One should notice that 



in these studies only Refs. |20j, |60j consider the recently reported / (1790). In a different 
fashion within the quark model, the a (980) and / (980) mesons were analyzed in Ref. f6ll | . 
being the effect of the two-pseudoscalar meson thresholds the responsible for the substantial 
shift to a lower mass than what is naively expected from the qq component alone. This gives 
rise to an important KK and nrj' components in the ao(980) and KK, rji], rj'r]' and rjf)' in 
the /o(980). 

The structure of the scalar mesons a (980) and /o(980) have been also investigated in 
the framework of a meson exchange based on the Julich potential model [H^ for im and nr] 
scattering. Whereas the / (980) appears to be a KK bound state the a (980) was found to 
be a dynamically generated threshold effect. Similar conclusions have been obtained in a 
chiral unitary coupled channel approach, where the / (600), the a (980), and the A"q(800) 
rise up as dynamically generated resonances, while the /o(980) is a combination of a strong 
S— wave meson-meson unitarity effect and a preexisting singlet resonance [H^. In Ref. 
|fiH van Beveren et al. describe the light scalar mesons as resonances and bound states 
characterized by complex singularities of the scattering amplitude. 

Finally, let us stress that Ref. [2(j presents an interpretation of the scalar mesons in a 
model constrained by the description of other hadron sectors. The same mechanism has 
been applied here to disentangle the structure of the new D s and D s j resonances. It drives 
to a final scenario that it is compatible with some other models in the literature and it differs 
from other results. Being the set of data so huge, and sometimes so poor, one always may 
find a positive or negative interpretation of some of them. Therefore, the final answer could 
only be obtained from precise experimental data that would allow to discriminate between 
the predictions of different theoretical models [65| . 



IV. SUMMARY 

As a summary, we have obtained a rather satisfactory description of the positive par- 
ity open-charm mesons in terms of two- and four-quark configurations. The mixing be- 
tween these two components is responsible for the unexpected low mass and widths of the 
D*j(2317), _Dsj(2460), and _Dq(2308). The same mechanism has been used to account for 
the spectroscopic properties of the light-scalar mesons. The obtained electromagnetic decay 
widths give hints that would help in distinguishing the nature of these states. We predict 
a ratio D s j(2A60) — > D^ , y/D s j(2A60) — > D* +/ y much larger than the one obtained in a pure 
qq scheme. We did not find any theoretical partner for the recently measured _Dj}(2632) 
whose existence awaits confirmation jsj. We encourage experimentalists on two different 
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directions: the measurement of the electromagnetic decay widths of the _D*j(2317) and the 
_D s j (2460), and the confirmation or refutation of the _Dj}(2632) that would help to clarify 
the exciting situation of the open-charm mesons. 

V. ACKNOWLEDGMENTS 

This work has been partially funded by Ministerio de Ciencia y Tecnologia under Contract 
No. FPA2004-05616, by Junta de Castilla y Leon under Contract No. SA-104/04, and by 
Generalitat Valenciana under Contract No. GV05/276. 



[1] BABAR Collaboration, B. Aubert et ai, Phys. Rev. Lett. 90, 242001 (2003). 
[2] CLEO Collaboration, D. Besson et ai, Phys. Rev. D 68, 032002 (2003). 
[3] Belle Collaboration, Y. Mikani et ai, Phys. Rev. Lett. 92, 012002 (2004). 
[4] S. Godfrey and R. Kokoski, Phys. Rev. D 43, 1679 (1991). D. Ebert, V.O. Galkin, and R.N. 
Faustov, Phys. Rev. D 57, 5663 (1998). M. Di Pierro and E. Eichten, Phys. Rev. D 64, 114004 
(2001). W. Lucha and F.F. Schoberl, Mod. Phys. Lett A 18, 2837 (2003). S. Godfrey, J. Phys. 
Conf. Ser. 9, 59 (2005). 
[5] Belle Collaboration, K. Abe et ai, Phys. Rev. D 69, 112002 (2004). 
[6] FOCUS Collaboration, J.M. Link et al, Phys. Lett. B 586, 11 (2004). 
[7] SELEX Collaboration, A.V. Evdokimov et ai, Phys. Rev. Lett. 93, 242001 (2004). 
[8] BABAR Collaboration, B. Aubert et ai, hep-ex/ 0408087, Belle Collaboration, B. Yabsley, 

|hep-ex/0507028| FOCUS Collaboration, R. Kutschke, E831-doc-701-v2. 
[9] E. van Beveren and G. Rupp, Phys. Rev. Lett. 91. 012003 (2003). 
[10] Yu.A. Simonov and J.A. Tjon, Phys. Rev. D 70, 114013 (2004). 
[11] R.N. Cahn and J.D. Jackson, Phys. Rev. D 68, 037502 (2003). 
[12] W.A. Bardeen, E.J. Eichten, and C.T. Hill, Phys. Rev. D 68, 054024 (2003). 
[13] D.S. Hwang and D.W. Kim, Phys. Lett. B 601, 137 (2004). 
[14] T. Barnes, FE. Close, and H.J. Lipkin, Phys. Rev. D 68, 054006 (2003). 
[15] A.P. Szczepaniak, Phys. Lett. B 567, 23 (2003). K. Terasaki, Phys. Rev. D 68, 011501(R) 
(2003). 

[16] HY. Cheng and W.S. Hou, Phys. Lett. B 566, 193 (2003). 

[17] Y.-Q. Chen and X.-Q. Li, Phys. Rev. Lett. 93, 232001 (2003). 

[18] T. Browder, S. Pakvasa, and A. Petrov, Phys. Lett. B 578, 365 (2004). 

[19] OS. Bali, Phys. Rev. D 68, 071501(R) (2003). 

[20] J. Vijande, A. Valcarce, F. Fernandez, and B. Silvestre-Brac, Phys. Rev. D 72, 034025 (2005). 
[21] J. Vijande, F. Fernandez, and A. Valcarce, J. Phys. G 31, 481 (2005). 
[22] Y. Suzuki and K. Varga, Lecture Notes in Physics M 54, 1 (1998). 

[23] J. Vijande, F. Fernandez, A. Valcarce, and B. Silvestre-Brac, Eur. Phys. J. A 19 383 (2004). 
[24] A. de Rvijula, H. Georgi, and S.L. Glashow, Phys. Rev. D 12, 147 (1975). 
[25] R.K. Bhaduri, LE. Cohler, and Y. Nogami, Phys. Rev. Lett. 44, 1369 (1980). 
[26] J. Weinstein and N. Isgur, Phys. Rev. D 27, 588 (1983). 

[27] S. Titard and F.J. Yndurain, Phys. Rev. D 51, 6348 (1995). A.M. Badalian and B.L.G. Bakker, 
Phys. Rev. D 62, 094031 (2000). 



11 



[28] F. Halzen, C. Olson, M.G. Olsson and M.L. Stong, Phys. Rev. D 47, 3013 (1993). 
[29] G.S. Bali, Phys. Rep. 343, 1 (2001). 

[30] SESAM Collaboration, G.S. Bali, H. Neff, T. Diissel, T. Lippert, and K. Schilling, Phys. Rev. 
D 71, 114513 (2005). 

[31] PW. Stephenson, Nucl. Phys. B 550, 427 (1999). P. de Forcrand and O. Philipsen, Phys. 

Lett. B 475, 280 (2000). 
[32] J. Vijande, P. Gonzalez, H. Garcilazo, and A. Valcarce, Phys. Rev. D 69, 074019 (2004). 
[33] P. Gonzalez, A. Valcarce, H. Garcilazo, and J. Vijande, Phys. Rev. D 68, 034007 (2003). 
[34] M.M. Brisudova, L. Burakovsky, and T. Goldman, Phys. Rev. D 61, 054013 (2000). 
[35] S. Eidelman et al, Phys. Lett. B 592, 1 (2004). 
[36] J. Hein et al, Phys. Rev. D 62, 074503 (2000). 

[37] UKQCD Collaboration, P. Boyle, Nucl. Phys. B (Proc. Supp.) 63, 314 (1998); et al, Nucl. 

Phys. B (Proc. Supp.) 53, 398 (1997). 
[38] UKQCD Collaboration, A. Dougall, et al Phys. Lett. B 569, 41 (2003). 
[39] N. Isgur and M. B. Wise, Phys. Lett. B 237, 527 (1990); ibid Phys. Lett. B 232, 113 (1989). 
[40] S. Godfrey, Phys. Lett. B 568, 254 (2003). 
[41] M. Nielsen, |hep-ph/0510277| 

[42] P. Colangelo and F. de Fazio, Phys. Lett. B 570, 180 (2003). P. Colangelo and F. de Fazio, 

and R. Ferrandes, Mod. Phys. Lett. A 19, 2083 (2004). 
[43] C. Chang, C.S. Kim, and G. Wang, Phys. Lett. B 623, 218 (2005). 

[44] T. Barnes, F.E. Close, J.J. Dudek, S. Godfrey, and E.S. Swanson, Phys. Lett. B 600, 223 
(2004). 

[45] F. E. Close and E. S. Swanson, Phys. Rev. D 72, 094004 (2005). 

[46] T. Mehen and R.P Springer, Phys. Rev. D 70, 074014 (2000). 

[47] P. Colangelo, F. De Fazio, and A. Ozpineci, Phys. Rev. D 72, 074004 (2005). 

[48] Belle Collaboration, K. Abe et al, |hep-ex/0507019[ 

[49] Belle Collaboration, S.-K. Choi et al, Phys. Rev. Lett. 94, 182002 (2005). 
[50] Belle Collaboration, K. Abe et al, |hep-e x/0512035 

[51] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 95, 142001 (2005). 

[52] CDF Collaboration, M.D. Corcoran, |hep~- ex/0506061 . 

[53] BES Collaboration, M. Ablikim, et al, Phys. Lett. B 607, 243 (2005). 

[54] V.V. Anisovich and A.V. Sarantsev, Eur. Phys. J. A 16, 229 (2003). 

[55] R.L. Jaffe, Phys. Rev. D 15, 267 (1977). 

[56] C.J. Morningstar and M. Peardon, Phys. Rev. D 60, 034509 (1999). 

[57] J. Weinstein and N. Isgur, Phys. Rev. Lett. 48, 659 (1982); ibidPhys. Rev. D 41, 2236 (1990). 
[58] C. McNeile and C. Michael, Phys. Rev. D 63, 114503 (2001); W. Lee and D. Weingarten, 

Phys. Rev. D 61, 014015 (2000); D.-M. Li, H. Yu, and Q.-X. Shen, Commun. Theor. Phys. 

34, 507 (2000); L.S. Celenza et al, Phys. Rev. C 61, 035201 (2000). 
[59] C. Amsler, Phys. Lett. B 541, 22 (2002); C. Amsler and F.E. Close, Phys. Rev. D 53, 295 

(1996); F.E. Close and A. Kirk, Eur. Phys. J. C 21, 531 (2001). 
[60] F.E. Close and Q. Zhao, Phys. Rev. D 71, 094022 (2005). 
[61] N.A. Tornqvist, Z. Phys. C 68, 647 (1995). 

[62] G. Janssen, B.C. Pearce, K. Holinde and J. Speth, Phys. Rev. D 52, 2690 (1995). 

[63] J. A. Oiler and E. Oset, Phys. Rev. D 60, 074023 (1999). J. A. Oiler, E. Oset and J.R. Pelaez, 

Phys. Rev. Lett. 80, 3452 (1998). 
[64] E. van Beveren, et al, Z. Phys. C 30, 615 (1986); E. van Beveren and G. Rupp, Eur. Phys. 



12 



J. C 10, 469 (1999). 

[65] C. Amsler and N.A. Tornqvist, Phys. Rep. 389, 61 (2004) and references therein. 



13 



TABLE I: cs and cn masses (QM), in MeV. Experimental data (Exp.) are taken from Ref. |3 
except for the state denoted by a dagger that has been taken from Ref. 0. 



nL J p 


State 


QM (cs) 


Exp. 


State 


QM (cn) 


Exp. 


IS 0" 


D s 


1981 


1968.5±0.6 


D 


1883 


1867.7±0.5 


15 1- 


D* s 


2112 


2112.4±0.7 


D* 


2010 


2008.9±0.5 


1P 0+ 


D* sJ (2317) 


2489 


2317.4±0.9 


L>S(2308) 


2465 


2308±17±15±28t 


IP 1+ 


D s j (2460) 


2578 


2459.3±1.3 


L>i(2420) 


2450 


2422. 2±1. 8 


IP 1+ 


D s i(2536) 


2543 


2535.3±0.6 


L>?(2430) 


2546 


2427 ± 26 ± 25 


IP 2+ 


D s2 (2573) 


2582 


2572.4±1.5 


D|(2460) 


2496 


2459±4 



TABLE II: cnsn and cnnn masses, in MeV. 



cnsn 


cnnn 


J p = 0+ 
1=01=1 


J p = 1+ 
1=0 1=1 


J p = 0+ 
I = 1/2 


2731 2699 


2841 2793 


2505 



TABLE III: Probabilities (P), in %, of the wave function components and masses (QM), in MeV, 
of the open-charm mesons once the mixing between qq and qqqq configurations is considered. 
Experimental data (Exp.) are taken from Ref. jiH] except for the state denoted by a dagger that 
has been taken from Ref. lal. 



7 = 


I = 1/2 


J p = 0+ 


J p = 1+ 


J p = 0+ 


QM 
Exp. 


2339 2847 
2317.4±0.9 - 


QM 
Exp. 


2421 2555 
2459. 3±1. 3 2535.3±0.6 


QM 
Exp. 


2241 2713 
2308±17±15±28 t - 


P(cnsn) 
P(cs l3p ) 
P(cs 2 3 P ) 


28 55 
71 25 
~ 1 20 


P(cnsn) 
P(cs 1 i P ) 
P(cs l3p ) 


25 ~ 1 
74 ~ 1 
~ 1 98 


P(cnnn) 
P(cnip) 
P(cn 2 p) 


46 49 
53 46 
~ 1 5 
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TABLE IV: Comparison of 90% C.L. limits on radiative transitions obtained by CLEO 2] and 
Belle 0| with our results (QM) and those of two different quark models, Refs. 0, l4o| . based 
only on qq components. The BR's are with respect to the decay D*j(2317) — > Dftr for the 
D*j(2317) and with respect to the decay Z) s j(2460) -> D* + ir° for the £> sJ (2460). We have assumed 
T(D* sJ (2317) -> D+tt°) « T(D sJ (2460) -» D* + ir°) w 10 keV as explained in the text. We have 
also quoted results obtained by light-cone QCD sum rules 47 1 and vector meson dominance [i^ . 



Transition 


Quark models 
QM Ref. [12] Ref. [40] 


Experiments 
CLEO [2] Belle [3] 


Other approaches 
Ref. [47] Ref. [42] 


^(2317) -» L»*+ 7 
^(2317) - £>+7 
D sJ {2AW) - D s *+7 
D s<7 (2460) -» D+7 


0.16 0.17 0.19 
0.0 0.0 0.0 
0.006 0.47 0.55 
0.67 0.51 0.62 


< 0.059 < 0.18 

< 0.052 < 0.05 

< 0.16 < 0.31 

< 0.49 0.55±0.13±0.08 


0.4 - 0.6 0.085 
0.0 0.0 
0.06-0.11 0.15 
1.9-2.9 0.33 



TABLE V: Masses (QM), in MeV, of the recently measured charmonium and B c states obtained 
within the model of Ref. I2lj| used in this work. 



Name Mass Ref. 


n 2S+1 Lj QM 


X(3940) 3943 ± 6 ± 6 [48] 

Y(3940) 3943 ±11 ±13 [49J 
A^ 2 (3940) 3931 ±4± 2 [50J 
Y(4260) 4260 ± 8 ± 2 [51J 
£ c (6287) 6287 ±4.8 ±1.1 [52] 


2 1 P 1 3923 
2 3 P 3878 
2 3 Pi 3915 
2 3 P 2 3936 
4 3 Si 4307 
l 1 ^ 6277 
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